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After a rapid increase in temperature across the deconfinement temperature T^, 

pure gauge theories exhibit unstable long wavelength fluctuations in the approach to 

. equilibrium. This phenomenon is analogous to spinodal decomposition observed in 
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■rj" . condensed matter physics, and also seen in models of disordered chiral condensate 
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o 



X 



formation. At high temperature, the unstable modes occur only in the range < A; 



occurs via spinodal decomposition for SU (2) at temperatures T > and for SU{3) 

_ . for T S> Td- For SU (3) at temperatures T > T^, nucleation may replace spinodal 

D . 

r~| , decomposition as the dominant equilibration mechanism. Monte Carlo simulations 



of SU{2) lattice gauge theory exhibit the predicted phenomena. The observed value 
of kc is in reasonable agreement with a value predicted from previous lattice mea- 
surements of mo- 
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I. INTRODUCTION 



In a heavy ion collision, changes in local energy densities of 1 GeV/ fm? and higher 
occur on a time scale less than 1 fm/c. During the initial stages of the formation of a 
quark-gluon plasma, the system is thermodynamically unstable. Some consequences of this 
instability can be studied with both simulation and analytical techniques. In particular, 
we find the existence of exponentially growing long wavelength modes in the approach to 
equilibrium. Such behavior is referred to as spindodal decomposition in condensed matter 
physics. [|T| [0 We will study spinodal decomposition in finite temperature gauge theories 
using both analytical and simulation techniques. The existence of exponentially growing 
long wavelength modes cutoff at a fixed wavelength depends only on the features of the 
equilibrium effective action in an unstable region. Although the methods we use are not 
based on the true real-time evolution of the system, they give a consistent picture of some 
basic features of real-time behavior. 

We consider first the generic case of a pure SU{N) gauge theory in which the temperature 
is raised rapidly from a temperature less than T^, the deconfinement temperature, to a 
temperature above T^. We refer to such a rapid increase in system temperature as a quench; 
this nomenclature is borrowed from statistical mechanics, where it is usually applied to rapid 
cooling of a system below a critical temperature. It is appropriate here because it is the 
high temperature phase of the pure SU{N) gauge theory which spontaneously breaks global 
Z{N) invariance. At temperatures T < T^, the Z{N) global symmetry associated with 
confinement is unbroken. The standard equilibrium order parameter for Z{N) symmetry 
breaking is the Polyakov loop L, defined in equilibrium as the path-ordered exponential 



The Polyakov loop can be defined for a general density matrix by its association with the 
projection operator onto gauge invariant states. Because of Z(A^) symmetry, {Lp) = 
jj {Trp^L)) vanishes below T^. When the temperature is rapidly increased to T > T^, 
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{Lp) =0 is no longer the stable state of the system, and must evolve to a new equilibrium 
state with {Lp) 7^ 0, reflecting the transition to the gluon plasma phase. We will use the 
Polyakov loop as the primary tool for the study of this process. 



II. EFFECTIVE POTENTIAL 

The instability of {Lp) =0 at high temperatures follows from the one-loop finite tem- 
perature effective potential for the Polyakov loop, as derived by Gross, Pisarski and Yaffe 
and Weiss. Q For simplicity, consider the case of SU{2). At any spacetime point, the 
Polyakov loop can be diagonalized, and we can write 



(2.1) 



where < g < 1. It is convenient to introduce the parameter which is related to q by 
q= (1 — '?/') /2. The order parameter is given by 

= Irri^L = cos[7r(l -^)/2] . (2.2) 

The effective potential at one loop for gauge bosons in a constant Polyakov loop background 
is 

nw = -^ + ^(i-^-f . (2.3) 

The one loop result dominates the effective potential for T ^ ^qcd due to asymptotic 
freedom. The first term is the standard black body result, obtained when = 1. Figure 1 
shows \^ as a function of if), normalized such that V{il) = 1) = 0. The use of the variable 
il) makes the Z{2) symmetry of the potential under ip —^/^ manifest. Note that the 
equilibrium value of ip is ±1, corresponding to Lp = ±1; = 0, corresponding to Lp = 0, 
is a maximum of V{iIj). 

Our picture of the quenching process is that the system is initially in a state where ip is 
equal to zero at some temperature below T^. When the system is quickly raised to a new 



temperature T > Td, the system is still in the state with = 0. However, the system 
is unstable, and must eventually find its way to either ip = +lor ip = —1. Because we 
quench into a region of the phase diagram where V"{tp) < 0, the system will decay to the 
equilibrium state via spinodal decomposition. 



III. LANGEVIN MODEL 



In order to study the dynamics of this transition, we use the effective action of Bhat- 
tacharya combined with Langevin dynamics. H The effective action takes the form 



The equilibrium distribution will be 



2r 



(3.1) 



exp [-Seff [i)]] . 



(3.2) 



We postulate Langevin dynamics of the form 



dip{x, t) 



-V- 
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where the white noise rj is normalized to 



+ ri{x,t) 



(3.3) 



{ri{x, r)7]{x\ t')) = 2T5%x - a;')5(r - r'). 



(3.4) 



Assuming translation- invariant initial conditions, the expectation value "^{t) = {ipix, r)) 
will evolve away from the unstable value \1' = approximately as 



>.rp3 



dr 3 
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(3.5) 



This gives initial exponential growth, which slows down as equilibrium is approached. 

The late-time relaxational behavior of ip is controlled by the Debye screening mass m^j, 
given by m|) = 2(yf^T^/3. Near equilibrium, any effects of initial conditions decay exponen- 



tially as exp 



27r2 rr / 7,2 



(k + mi))T . However, for early times, the initial conditions contribute 



to (^(k,T)ilj{—k,T)) a term of the form 



ip{k, 0)ij{-k,0)exp 



9' 



(3.6) 



where k1= = m\yj1. Modes with k < kc are initially not damped but grow exponen- 

tially, with the k = mode growing the fastest. This is a characteristic feature of spinodal 
decomposition with a non-conserved order parameter. 

Because pure SU{2) gauge theory has a second-order deconfining transition, spinodal 
decomposition will occur after quenching to any temperature T > T^. The situation is more 
subtle for SU{3). The one-loop effective potential is unstable at L^t^ = 0. For temperatures 
sufficiently large that the one-loop effective potential for L is reliable, spinodal decomposition 
will occur. However, the first-order character of the pure SU{3) transition implies the 
existence of a metastable phase with (Lp) = Ofor some range of temperature T ^ T^. This 
in turn implies that nucleation followed by bubble growth is the important mechanism for 
attaining equilibrium for some range of temperatures just above T^. 



IV. DYNAMICAL QUARKS 

When dynamical particles in the fundamental representation, e.g, quarks, are included, 
the Z{N) symmetry is explicitly broken. At low temperatures, Lp{x) ^ 0. However, after a 
rapid quench to T > T^, Lp must still change to its new equilibrium value. This change may 
involve spinodal decomposition or nucleation as well as relaxational processes. The relevant 
potential for SU{2) with Nj massless fermion flavors has the form 

= ^ (l - V^f + (7 + - i^') (1 - i^? . (4.1) 

With Nf = 2,V {tp) is negative for t/; < |, as shown in figure 2, so spinodal decomposition 
will take place for a significant range of initial conditions. |^ |^ Thus the initial value of 
Lp in the confined phase will determine whether spinodal decomposition occurs. These 
considerations are independent of any static critical behavior, e.g., the existence of a decon- 
finement or chiral phase transition for particular values of quark parameters. Because Z{N) 
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symmetry breaking terms in the effective potential vanish in the hmit of infinite quark mass, 
spinodal decomposition is to be expected for sufficiently large quark masses. 

V. MONTE CARLO RESULTS FOR SU{2) 

We have tested these theoretical ideas with simulations of rapidly quenched pure SU{2) 
gauge theory. Lattices of size 32^ x 4 and 64^ x 4 were equilibrated at f3 = 2.0; the decon- 
finement transition for Nt = 4 occurs at (3d = 2.2986 ± 0.0006. [0 The coupling constant 
was increased instantaneously to (3 = 3.0, and the approach to equilibrium monitored via 
the Polyakov loop and other observables. For simplicity, the heat bath algorithm was used 
both for equilibration at low temperature before the quench and for subsequent dynamical 
evolution after the quench. While this time evolution is of course not the true time evolu- 
tion of the non-equilibrium quantum field theory, features such as spinodal decomposition 
which depend only on the equilibrium effective action will occur with any local updating 
algorithm which converges to the equilibrium distribution. The abrupt change in /3 is a po- 
tential cause of concern with this procedure, since the lattice spacing, and hence the physical 
volume, changes in all directions when (3 is changed. However, the large spatial sizes used 
should mitigate this effect. It is also possible to maintain a constant spatial volume with 
anisotropic lattice couplings. As an alternative, we have also taken 32^ x 4 sections from 
32^ X 10 configurations equilibrated at /? = 3.0, with similar results. However, the initial 
configurations produced in this way suffer from an obvious problem due to abrupt joining 
of the new boundaries in the time direction. 

In figure 3, we show the Polyakov loop expectation value versus Monte Carlo time for 
one 64^ x 4 simulation. Comparison with a numerical integration of equation 3.5 shows 
qualitative agreement, although the equilibrium value of (Lp) measured in simulations has 
a multiplicative renormalization. Figure 4 shows the Fourier transform of the connected 
Polyakov two point function S{k, r) for low values of the wave number as a function of Monte 
Carlo time for the same simulation. Note the early exponential rise in these modes, followed 
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by a sharp disappearance as the Polyakov loop reaches its equihbrium value, characteristic 
of spinodal decomposition. Only the low momentum modes exhibit this growth; above kc 
no such growth occurs. Although the general behavior of S{k, r) is the same for each run, 
many details are run dependent. In this particular run, the k/T = 0.68 mode achieves a 
larger amplitude than k/T = 0.56, which is atypical. In some runs, there is clear evidence 
for mode-mode coupling, reflecting the nonlinearity of the system. For each run, we have 
estimated the rate of growth of each low-momentum mode by fitting log{S{k, r)) to a straight 
line in r for early times. We can extract kc as the value where the growth rate is zero. From 
equation 3.6, the growth rate of each line is proportional to k1 — k"^ for the linearized theory, 
but this may not represent the true time evolution. In any case, the growth rates measured 
in each run are highly dependent on initial conditions. In figure 5, we plot these growth rates 
versus /c^/T^ for the same run used in figure 4. The error bars are naive estimates of the 
error for each growth rate for this particular run. The x intercept provides an estimate of k^ 
for each run. Using multiple 64^ runs aX (3 = 3.0, we have estimated k^/T to be 1.09 ± 0.08. 
The principle errors in this estimate come from the sensitivity of individual modes to initial 
conditions and the discrete character of k on the lattice. We can compare this with lattice 
measurements of the Debye screening length, assuming the one-loop relation 

kc tTioiT) 

= 7- 5.1 

holds in general. Using the results of Heller et al. for m£,{T), we obtain kc/T = 1.35(5). 
We consider this to be reasonable agreement, given the many uncertainties involved. 

Somewhat less sensitive to initial conditions is the traditional observable k*, which is 
defined to be 

, */ N _ /o°° dk k S{k, t) 

" ~ C dk S{k, t) ■ ^^-^^ 

It is convenient to plot T/k*{T) versus r, as in figure 6. Its early time behavior is predicted 
from equation 3.6 to be 1/A;*(r) ~ ^/T, typical of length scales in diffusive processes. We 
do observe this behavior, but only for very early times; on 64^ x 4 lattices, this occurs for 



< r ^ 120. For the run shown in figure 6, T/k*{T) is fit well by ar^ for < r < 120, 
where p = 0.54 ± 0.03. 

VI. CONCLUSIONS 

As we have shown above, theoretical and simulation results both indicate the relevance 
of spinodal decomposition in the equilibration of a gluon plasma after a rapid quench to 
high temperature. Since simulations and theoretical analysis were both carried out in Eu- 
clidean space, the significance of these results for experiment is not clear. However, spinodal 
decomposition is a general phenomenon whenever a system is in an unstable initial state, 
as determined by some local variable. Since we have considered a rapid heating of a finite 
temperature gauge theory, it is natural to ask what happens when rapid cooling takes place, 
as might occur in the late stages of the expansion of a quark-gluon plasma, or in the early 
universe. At low temperatures, we expect that Vltp) has a single minimum and V"{iIj) > 
everywhere; this is required by Z[N) symmetry in the case of a pure gauge theory. Thus the- 
ory predicts the absence of unstable modes, and that relaxational processes should dominate 
the approach to equilibrium. We have performed simulations of such a cooling process, in 
which 32^ X 4 lattice configurations equilibrated at /? = 3.0 are suddenly cooled to P = 2.0. 
Examination of the data from these runs shows no sign of spinodal decomposition. 

It is interesting to contrast the phenomenon of spinodal decomposition in a rapid heating 
with the formation of a disordered chiral condensate (DCC) in a rapid cooling. ||^ In both 
cases, exponentially growing low-momentum modes occur when an unstable initial state 
must equlibrate to a final state associated with a broken symmetry. In the case of a DCC, 
chiral symmetry is the relevant symmetry. For massless quarks, the symmetry is sponta- 
neously broken at low temperatures. With sufficiently light quarks, the overall structure of 
the potential survives, and the approach to equilibrium is fundamentally the same as the case 
of massless quarks. In the absence of fundamental representation particles such as quarks, 
Z{N) symmetry breaks spontaneously at high temperatures. Provided there are sufficiently 
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few light quarks, spinodal decomposition still occurs at sufficiently high temperature. 

Just as the formation of a DCC may lead to enhanced production of low-momentum pi- 
ous, spinodal decomposition should lead to enhanced production of low-momentum gluons 
in the early stages of plasma formation. The characteristic scale for such a phenomenon 
would be gT. As the system equilibrates, this enhancement of the small-A; part of the gluon 
distribution will be eliminated. It is interesting to note that another, perhaps related, mech- 
anism for the creation of exponentially growing low-momentum modes has been proposed; 
it has been argued that such instabilities can significantly affect the temporal evolution of 



the system. [10 



In contrast to the SU{2) results presented here, the case of SU{3) presents additional 
complications due to the presence of the confined phase as a metastable state over some range 
of temperature in the deconfining regime. Study of the SU (3) dynamics may be helpful in 
exploring the limits of met ast ability. It would also be of interest to study via simulation the 
effect of unquenched quarks on the dynamics. Because the dynamic part of the simulation 
is much faster than creating equilibrated field configurations, this could be done fairly easily 
if large equilibrated unquenched lattice field configurations became available. 
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FIGURES 



FIG. 1. The effective potential F(V')/r^ for SU{2) with ATj = as a function of tp. 
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FIG. 2. The effective potential F(V')/r^ for SU{2) with AT/ = 2 as a function of 
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FIG. 3. Polyakov loop as a function of Monte Carlo time. 
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FIG. 4. S{k,T) versus Monte Carlo time. 
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FIG. 5. Growth rate of S{k, r) versus k'^/T'^. 
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FIG. 6. T/k*{T) versus Monte Carlo time. 
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